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NONLIFTING WING-BODY COMBINATIONS WITH CERTAIN GEOMETRIC RESTRAINTS HAVING 

MINIMUM WAVE DRAG AT LOW SUPERSONIC SPEEDS 1 


By Har vard Losiax 


SUMMARY 

Several variational problems involving optimum wing and 
body combinations having minimum wave drags for different 
kinds of geometrical restraints are analyzed. Particular 
attention is paid to the effect on the wave drag of shortening the 
fuselage and, for slender axially symmetric bodies, the effect 
of fixing the fuselage diameter at several points or even of fixing 
whole portions of its shape. 

INTRODUCTION 

Recently several authors have used linearized theory to 
study the wave drag of wing-body combinations traveling 
at supersonic speeds (see, o. g., refs. 1 to 5). These studies 
have clearly demonstrated the importance of finding the 
wave drag of a whole airplane rather than the separate 
wave drags of its various parts (wings, fuselages, etc.), 
since the magnitude of the interference terms can pre- 
dominate. In effect, this means that various optimization 
problems for bodies — such as the problem of finding the body 
shape having a minimum wave drag for a given volume — 
should be re-examined when interfering wings or other bodies 
are in the same flow field. In many cases the solution to 
the new problem differs from the body-alone problem only 
in interpretation. 

The purpose of this report is to study minimum wave-drag 
combinations which satisfy a few of the many possible 
geometric restraints pertinent to the interests of airplane 
designers. An attempt has been made to analyze the 
various problems in a unified manner so that extensions to 
other kinds of restraints can bo deduced. 


LIST OF IMPORTANT SYMBOLS 
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aspect ratio 

source distribution equivalent to wing in 
sense defined by equation (3) 
multipole distribution of order n 
wave drag 

portion of drag due to all the nth order 
multipoles for n f> 0 
See equation (8). 

additional drag resulting from restraint (See 
eq. (11).) 

restraints defined in equations (19) 
distance between apexes on x axis of forecone 
and aftercone enclosing wing (See fig. 3.) 
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length of basic body 

length of modification to basic body 

Mach number 

PJJJ 

2 

average body radius 

fuselage area in cross section normal to the 
free stream 

normal (to free stream) projection of wing 
area in section cut by plane x^xpPy, cos 6 
(See fig. 2.) 
speed of free stream 
volume 

Carteshnr coordinate system (See fig. 1.) 
source distribution representing the fuselage 
modifications 

polar coordinate (See fig. 1.) 
free-stream density 
See equation (17). 
velocity potential 

BASIC THEORY AND ASSUMPTIONS 


BASIC THEORY 

Many of the discussions and derivations contained in the 
following are carried out on the assumption that the reader 
is familiar with the concepts presented in reference 4 which 
should be considered as a first part, to this report. In 
particular, an acquaintance with the solutions to the wave 
equation referred to as “multipoles” is assumed, together 
with Hayes’ invariance principle and the consequent multi- 
pole distributions equivalent to a wing in the sense that 
both induce the same momentum spectrum at infinity. 

The entire analysis used herein is based on the assumptions 
and idealizations necessary to develop the linearized equation 
for the velocity potential, <p, in supersonic flow, namely 

fiVn— <Pn~ <P«=0 (1) 

where 0 2 =M 2 — 1 and the reference coordinate system 1 is 
shown in figure 1. The analysis is further restricted to the 
solution of problems involving a given uncambered and un- 
twisted wing mounted centrally on a vertically symmetrical 
fuselage, the entire configuration being at zero angle of 
attack. 


1 Supersedes NACA TN 3C67, “Wing-Body Combinations With Certain Geometrio Restraints Haring Low Zero-Lift Wave Drag at Low Supersonic Mach Xumbere,” by Harv ard Lomax 


lose. 

' It should b« stressed that the rails Is parallel to the free-stream direction (wind axes). 
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Figure 1. — Coordinate systems used in analysis. 

ADDITIONAL ASSUMPTIONS 

We now make the two additional assumptions: one, the 
value of pA, where A is the wing aspect ratio, is small; and 
two, the value of p EJL e , where R is the average body radius 
and 2 L, is the distance along the x axis in which the multipole 
strengths differ significantly from zero, is small. 

One can evaluate the significance of these assumptions by 
studying their implications relative to the multipole distri- 
butions used to simulate the "wing and body. Suppose, for 
example, a group of nth order multipoles is placed along the 
body center line, their strengths,' C n ( x), being fixed by the 
condition that a circular cylinder in the vicinity of the body is 
a stream surface when the velocity field induced by these 
multipoles is combined with the velocity field induced by the 
source sheets representing the given wing. With the assump- 
tions of small pA and pR/L,, mentioned above, the C„’ s, for n 
greater than 0, can be shown (see, e. g., ref. 4) to have a negli- 
gible effect on the wave drag. Hence, all the multipoles (for 
0) that combine with the wing to'make a circular cylinder 
a stream surface and any additional multipoles (for n>0) 
added to make the body have mild distortions from such a 
surface are negligible in evaluating the wave drag. There- 
fore, under the assumptions mentiond above, out of all the 
singularities distributed along the body axis, it is necessary^ 
in studying the wave drag, to consider only the sources 
(multipoles for which »=0). 

With the restrictions to small values of pRJL „ and mild 
body distortion (see Ward, ref. 6, for a discussion of orders of 
magnitude), slender-body theory can be used to calculate the 
body shape, and on the basis of this theory one can show 
(see ref. 4, Appendix B) that S/(x), the body cross-sectional 
area measured normal to the free stream, is completely deters 
mined by the axial source distribution alone. Hence, if only 


the exposed .panels of the wing are used to calculate the £7»’s, 

Co is negligible and the entire axial source distribution a 0 (x) ■ 
is related to the geometrical properties of the body by the 
relation 

a*(x)^U 0 ^=U a S/(x) ( 2 ) 

THE WING EQUIVALENT SOURCE DISTRIBUTION AND THE OPTIMUM 
CANCELLATION SOURCES 

Let the given wing lie in the Zt= 0 plane. According to 
Hayes’ theorem (ref. 7), the wing equivalent source distribu- 
tion [oo(a:)]tt is obtained by accumulating on the Xi axis, at 
a distance x from the origin, all the wing sources intercepted 
by the line x^x+Py&os 6, and then, for a fixed x, averaging 
these values as d varies from 0 to 2ir. Thus using thin-airfoil 
theory to relate the planar source sheet to wing geometry, 
one finds 

-U a«(x)=± J* SJ(x,0)dO (3) 

where S a ' (x,8)=b[dx[S a (x,9)] and S w (x,0) is the normal 
(to the x axis) projection of the wing cross-sectional area 
intercepted by the plane 3 a^r-f jS^cos 6 as shown in figure 
2. Without the addition of further restraints, the optimum 
source distribution along the x x axis is that which just cancels 
the wing equivalent source distribution. Further, this can 
be interpreted directly in terms of both fuselage and wing 


S(x,B ) = Normal projection 
of wing area along A A 



* The true oblique plane Is given by the equation 

zi^j-fftrtcos *+0naln 9 

but, to be consistent with the assumptions basic to linearised theory, the variation with i\ 
is neglected. 



WING-BODY COMBINATIONS WITH CERTAIN GEOMETRIC RESTRAINTS HAVING MINIMUM WAVE DRAG 115 


geometry by means of equations (2) and (3). Thus, with no 
further restraints, the best fuselage shaping, for a wing-body 
combination satisfying the assumptions discussed above, 
satisfies the equation 

S/(x)=~J* S a ' ( x,0)dO (4) 

and has any reasonably smooth cross-sectional form. Notice 
that the total volume taken out of the fuselage is exactly 
equal to the total volume of the exposed portion of the wing. 
Hence, the total volume of the modified combination is the 
same as that of the original smooth cylinder. 

THE DRAG 

The total wave drag of a system can be expressed in terms 
of its actual or equivalent multipole distributions as 

Z?=2Do+S®* (5) 

i 

where D n is the drag -caused by the nth order multipoles 
a n (x) and is given by the equation 

f=~4 SbS-W^S-W I 

for n=0,l,2 ... (6) 

whoroffi„ <n ' H) (x) represents (d/hr) n+1 « 2 «(r). Under the assump- 
tions given above, the magnitude of X) D n is small. Let us 
designate it by D lt so that, in general, 

D=2Do~\~I-), (7) 

On the other hand, the total wave drag of a system com- 
posed of the combination of a wing and a body can also be 
written symbolically as 

D=D a -{-2D a t,-\-Di J (8) 

where 

D w drag of the wing alone 
D 6 drag of the body alone 
2D wb interference drag 

The various components of wave drag defined in equations 
(7) and (8) help one to evaluate more readily the drag reduc- 
tions that can be realized from appropriate fuselage shapings. 
Thus, if the fuselage shape satisfies equation (4), the total 
wave drag of the combination under the assumptions that 
p A and pR/L t are small can be written either as 


D=D. 

(9) 

or as 


D—Dv, — Db 

(10) 

If, in finding the fuselage shape, 



(a) the multipoles representing a wing and a body flying 

separately are assumed to represent the same 
wing and body when combined (i. e., the shape 
fields can be superimposed) , 

(b) the multipoles representing the fuselage are equal 

in magnitude but opposite in sign to the wing 
equivalent multipoles, 


then equation (10) holds without the assumption of small 
pA and PR/L,. 

In subsequent problems we will discuss the effects on the 
wave drag and fuselage area distribution of adding certain 
additional restraints to the body geometry. The addition 
of such restraints may or may not change the relation 
given by equation (10), but they must always add a term to 
equation (9) so that 

D—D.A-Drt' | 

y (id 

D„> 0 J 

WINGS CENTRALLY MOUNTED ON SLENDER 
QUASICYLINDERS 

This section is devoted to the solution of two problems 
involving a given uncambered and, untwisted wing mounted 
centrally at zero angle of attack on a tube that is cylindrical 
forward of some point ahead of the wing. The problems are, 
in both cases, to find the area distribution of the fuselage 
behind the cylindrical portion that will minimize the wave 
drag of the combination. 

SHORTENING THE FUSELAGE 

Remembering the assumptions fisted at the beginning of 
this section, let us consider the following problem: 

(i) Given a wing and a slender fuselage having the 
same normal cross-sectional area in all planes ahead of 
the plane x~ —I*/ (see fig. 3), what is the optimum fuse- 
lage area distribution behind the plane x=—Lo' if the 
fuselage must end at the plane x=U ? 

Of course, if lo>Lo (i. e., the body modification can extend 
over the entire range enclosed by the forecone and aftercone 
enclosing the wing) , the solution is already given by equation 
(4). Hence, in the following, Z 0 <C-£t>- 
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For simplicity of notations, let <^o(x) represent the sources 
along the fueslage center line and a 0 (x) represent the wing 
equivalent source distribution. Then, according to equation 
( 6 ) 


Do 1 



j*° w d^[o 0 , (r 1 )+ 


V(*i)][a 0 , (2^+Vte)]l^Zi-^l (12) 

where from the conditions stated in the problem and the 
geometric interpretation to the fuselage sources given by 
equation (2), 00 ( 2 ) is zero for values of x outside the interval 4 

-Lo'<x<k. 

Consider now a variation of equation (12) for a fixed 
a 0 (r) in the interval — L/ <x<Lo and a free variation of 
«o(x) in the subinterval — 7V<«<4- There results 

] 9 ( r* r* jj 

- 5D ° = ~ 4 ^ j ~ 2 J-w dxiWM [_J_£ 0 - ao'{x 2 )ln\x 1 —x^d^+ 


Using this expression, one can show that equations (8) and 
(12) yield 


Dub 

2 


b if'xWiif,,® 

U cd %) ^0 */ — ^*0 " 2*2 

U cd J “-^0 J ~W 


so that, by equation (13) 
D w & 1 


2?r(7„ 


J-W J-W Xi—z? 


Db 

2 


Hence, for any combination satisfying equation (14), once 
again . 

D=D K — Dh 

On the other hand, D r „ the increase in drag caused by 
shortening the fuselage, can also be obtained. Integrating 
equation (12) by parts, one has (note D„—2D 0 ) 


a*|<Wj=0 


D T , 1 
2 2 -kUJ 


taofe) +«o(*0] dx , 


X 1 —Z 2 


Integrate once by parts with respect to (since the varia- 
tions 5ao( — Lq) and 5ao(/ 0 ) must be zero). Then, by the 
fundamental lemma of the calculus of variations, the brack- 
eted term must be zero for — Zo'^x^Zo and one finds the 
condition 



do '(xjdx 3 |% 

Xi—x? + J- V 


\ 


tt o / (x 2 )cfo 3 _ 0 


— W ^Xi<l 0 (13) 


Equation (13) is an integral equation which can be inverted 
(by methods such as those outlined in ref. 8). Inverting, 
integrating once with respect to x, and expressing ao(x) 
and do ( 2 ) by means of equations (2) and (3), respectively, 
one finds 



S w '(x,0)d0-[ 


V(Zo — x)(V+x) 

2-r 2 



S a '(x u 9) 

(x t — x)V(V+ari)(x,— Z 0 ) 


(14) 


Combined with equation (13), this becomes 
D r , 




2 2v77 ra 

The derivative of equation (14) with respect to x gives 


Oo'(z)+“o'(a;)= 


— . r 

0 — x)(Lq' -{-x)J l ( 


x ° ao / (xi)V(L 0 , +Xi)(x 1 — Z 0 j 


wui 

—Lo'<x<l 0 


so 


A, 

2 


~2rUjft 0 


X— Xj 


fe)dxa 


Xj— Xj 


1 Ho 

’rJ-'V(2 1 

J, 


(Zx, 


x?)^(Lo -\~x%)(lo — X 3 ) 


L 0 ao / (x 3 )V(Lo , +X3)(x3— To) 


Xj— x 3 


(Zxj 


■] 


cZxj 


which reduces to 


which gives the optimum fuselage area distribution under the 
conditions and assumptions posed. 

The wave drag of the combination represented by equation 
(14) can be expressed either in the terms defined in equation 
(8) or (11). Let us first consider the form given by equation 
(8). If the expression for the drag of an nth order multipole 
distribution is integrated once by parts, there results since 
n^» +l) (-Lo')=a, ( » +1 >(Z o )=0 


D* 

2 


-&-S-w a -'” {xMx, SX 


On 0 ' +1} (Xj)dX3 
Xi— Xj 


D t , 


2 2irU, 


or, alternately, 

Dr 
2 


(16a) 


~2 Jj m *h L 


In 


|xi x^KZq+Lq') 


[V(Xi+40(X2-Zo) + V^+io'K*!— (o)] 2 


dxidxj (16b) 


CONSTRAINED FUSELAGE AREAS 


4 It Is necessary, for equation (6) to be valid, that LJ) and ^(Za) be rero. This implies 

that ao(r) must be continuous and if the body shape is given by equation (2), this, in turn 
means that the streamwise gradient of body cross-sectional area must be continuous. It was 
pointed out in reference 4 that oo(— Lf) and <zo(£q) will both be zero if the wing has no blunt 
edges along which thje normal component of the freestream Mach number is unity or greater. 


Another class of problems is that in which the magnitude 
of the fuselage area is fixed at various points. Suppose, for 
example, that a fuselage shaped according to equation (4) 
had in some region a cross-sectional area too small to be 



WING-BODY COMBINATIONS WITH CERTAIN GEOMETRIC RESTRAINTS HAYING MINIMUM WAVE DRAG 


117 



Figure 4. — Wing mounted on fuselage restrained at three sections. 

acceptable for some practical purpose. The question is, 
then, what is the best shape for given values of minimum 
fuselage cross-section area at given planes and what is the 
penalty in wave drag caused by the added constraints? 
Before considering the general case of an arbitrary number 
of restraints, let us first consider the simple problem: 

(ii) Given a wing, what (under the various assump- 
. tions given above) is the area distribution of the ad- 
joining fuselage which has a prescribed area at three 
given stations (the initial, the final, and an intermediate - 
station x—dy, see fig. 4) and yields a minimum wave 
drag for the combination? 

As before, let Oq(x) represent the wing equivalent source 
distribution. Then the drag caused by the restraints can 
bo written 




ao'(&)+U IB S/'(x a ) 


Xi—Za 


dx 2 
(16) 


where S f (x ) is the unknown fuselage area to be optimized. 
For simplicity, replace the unknown S r (x) by a(x) where 


Lot 


j(x)= J-L 0 IT d ^+ s ^ == t f s^de+SAti 

<r( — Lq)=So 
<t(Lo)=S 2 

a- 0 (x)=<r(x); ~L Q <z<di 
a 1 (x)=cr(x); di<x<Lo 


(17) 


(18) 


and the restraints on the fuselage area give the relations 

<r 0 , (z)cfo=[Si-6'o+^ iSUdi, »)<&]= Jo (19a> 

<r 1 / (r)dr=[~iS , 2 -S I -^ : (19b) 

where Jo and Jy are constants fixed by the given con- 
straints. Notice 

J^Si-So) (19c) 


so the constant Ji is a measure of the difference between 
the initial and final areas. 

Using the usual variational techniques, we can write, for 
the quantity to be minimized, 


Drt . Pa^O 

q 4-s - 




or 




By taking the variation and satisfying the conditions at the 
end points, one obtains the two simultaneous integral equa- 
tions 


C d 1 g-p "(xJdXj ■ C L » «y"{xjdxi _ V 
J-Xo Xy—X* ’’’Jtf, Xy—X* 2’ 

r d i <r 0 ' r (x 2 )dx 2 r r '» <Ty"(x i )ckh _ 

J-Z 0 Xy — Xi Jtf, Xy — Xi 2’ 


— L 0 <a: 1 <d, 
dy<^Xy<^La 


y ( 2 o> 


The set of equations (20) is identical to that analyzed by 
Adams (ref. 9) for bodies of revolution with fixed areas at the 
initial, final, and an intermediate section. In the interest 
of subsequent generalization, however, we will consider its 
solution in the following way: First write the equations (20) 
in the equivalent form 


cr "(x^dXj 

-L 0 Xy-X* 


One can show that 


7T> — Lo<^Xy<^dy 
dy<^Xy<^Lo 




A+Bx? 
V-L) S — X2 2 


Gy cosh 1 


Jp* — dyXy 

Lo\x* — di| 


( 21 ) 


( 22 ) 


is the solution to the integral equation (where A, B, and Cy 
are constants) since 


C L ° a”{x 2 )dx. 2 _ I 
J-Z 0 Xy—Xj | 


-.[u+ftoos-^)} 


— Lo<CXy <C.dy 
(23) 


dy<C.Xy<C.Lo 
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Since a'{ — Lo)=cr(Lo)=0 

A=£WW-di a 


which satisfies equation (21). The constraints can now be 
satisfied by means of the equations 


/(a:)=J ff n (x^dx l =A cos 1 ?+ 

C, [(*— af) cosh" 1 003-1 (if)] 


(24) 


and 


<r(x)—S 0 =A [x cos 1 + V-V — J 2 

V 003-1 (zf )] - ! 1 [} dl ~ x)2 cosh " 1 z^-di! - 


■yjLo 1 — d^idi— 2x) cos 1 (yj^) + V( Ai ! — di 2 ) (io 2 a^TJ 


(25) 


and 


Jo= — ^ [di V Lt?—di+Lt? cos " 1 (jX~)]+ 

Tdi cos -1 W-di s ] 


Ji=— ^ irdi^Ltf—di 2 


Solving for (7i and 5, we find 

B ~ ^w=d?r { 1,d ‘ J ^ J ' [* cos ~‘ (ir) + V£o '-*’]} 




in 2 COS -1 


G?)]} 


From equation (25) the fuselage cross-sectional area can he written 
S f (x)=So—^. j* T S a (xft)dB 1 ^ L J_ d 2 y {Vv-diz) V(£o*— **) (LJ-dd)- 


U'b-dtf cosh- 1 ^^]+ Jl {( i o 2_ ^ I)S c°s 1 (io1 ^ t^ 0 * * )+ 

(Lo s —dix) cos” 1 (^^-{-(x-d,) 2 [diV^ 2 -di 2 +io 2 cos" 1 cosh" 1 ^Z^|]} ) ( 2 °) 


In terms of the wing, body, and interference drag components defined in equation (8), the total wave drag is 


E = Q>J2±+(S 0 -SJ B+Ci [mSx-So cos- 1 cos 1 (-^)] 


(27a) 


2 2 2 

where B and Oi are given above. The equation for D„ is 


* j j : (^X 0 Vo 2 -2xJ 0 J, rd 1 Vio 2 -d 1 ! +io 2 COS 1 (^)]+ 
g 7r(iio — ctx J \ L. x J 

J 1 *{lJ-d 1 *+2d^fLJ=d? cos" 1 (^)+V [° 03-1 (^)J }) < 27b ) 


If the additional specification, is made that the initial and 
final areas are the same, the solution simplifies considerably, 
since, for such cases, Ji=0 and equations (26) and (27) 
reduce to 


SA)=s.-LJ”s.(*,e)i»+( Z7 ^ 3?? [s.-s«+ 

h£ r Sa(dl,e) [w-**) VCV-a^CV-di 2 )- 


(28) 




and 


D T 


2 (Zd 


(29b> 


Often the exact statement of the restraint is that S(x) 
shall not be lees than S t at x=di. In such cases care must 
be used in applying equations (26) and (27) or (28) and (29) , 
since they are only valid when the fuselage cross-sectional 
area at di is exactly If such is the case, equations (26) 
and (28) give the optimum body shape only if Jo>0, that is. 
only if 


& > J o * r 


(29a) 


Otherwise the opt imum variation of area is given by equation 

( 4 ). 

Next let us generalize the analysis leading to equation 
(26) and (27) by considering the following problem: 
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(iii) Given a wing, what is the area •distribution of 
the adjoining fuselage' which has prescribed areas at 
stations (including the initial tail-final ones fixed 
by the Mach forecone and aftercone enveloping the 
■wing, see fig. 5) and yields a minimum wave drag for 
the combination? 



and 

S w (x,6)dB - 1 

W-dr [d, cos-' j- (34) 

The wave drag due to the restraints can be obtained by 
using equations (32) and (16). Thus 

^=lBU 0 -a n+1 )+l± G t [w,-. 0 cos-'(|)- 

<L,+i cos -1 (35a) 

or in terms of the components defined in equation (8) 

f-y-T W s ‘- s -^+i( 0 ‘ [•*-& «•- d)~ 

S n+ i cos"' (^)] ’ (35b) 

where c t =a{di). Notice <r 0 =jS 0 and a n+ i—S n+1 , so when 

So~S n+ i, 

™=£±) Ct^-Sv) (35c) 

2 2 

.or 

t -So) (35d) 

2 2 2 f=i 

Finally, using the known values of S f ( x) at d„, v=0, 1, 
. . . , w+1, one obtains the n+1, simultaneous equations 


Fiqoee 5. — Wing mounted on fuselage restrained at n+2 sections. 


By analogy with equation (22), the integral equation 
for a"(x ) (where <r(a;) is defined by eq. (17)) that must be 
satisfied for a minimum wave drag can be written at once 
in the form 


J- 


u «>> 


S'jXg) dXy 


df\3u < \d(.{. 1,7- 0,1, • . *,7 v (30) 


l-r. aii— 

where d 0 =— Xo and d„ +1 =Zo. The quantity 


(31) 


is a solution to equation (30) since it yields 

Sl^-~i B +V °< (rr)-? « 


in which 


d,_!<a;j<d„ r=l,2, . . .,77+1 (32) 



Further, it is apparent from equations (24) and (25) that, 
with the conditions <r'(— L o )=a'(L o )=0 

cr'^—syz^+i) Ctidi-X) cosh- 1 (33) 

480104— 08 9 U ' 


(i5-i§o)=-^J o 2r S w (d v ,0)d0~R [c?^/Ly=d7+ 

V cos- 1 (^7)]-|S^{(^-d») s 

[* “ c "(x ! )+ 

y/V-d^Jj-; u=l,2, . . ., Ti+l (36) 

which determine the 77+ 1 constants B, Gi, Cs, . . ., C n . 
These, in turn, fix the shape, through equation (34), and 
the wave drag, through equations (35). 

Solutions similar to the above are presented in references 
10 and 11, and are used therein to calculate the drag of bodies 
of revolution having their areas specified at a given number 
of stations. Such a method has the advantage of giving the 
lower bound to the drag of bodies whose areas have been 
measured at a discrete number of places and, further, of 
giving a value representative of all area variations in the 
vicinity of the calculated optimum. Reference 10 contains 
a tabulation of the constants necessary to evaluate the 
minimum drag of an area distribution fixed at 19 points. 

WINGS CENTRALLY MOUNTED ON SLENDER CLOSED 
BODIES OF REVOLUTION 

In the preceding section the interference between the 
central portion of the airplane and its nose and tail regions 
was neglected. In this portion we will consider the entire 
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fuselage, assuming, first, it is a slender closed body and, 
second, it can be calculated in the presence of the wing, 
using the same postulates given in the previous section 
“Basic Theory and Assumptions.” 

UNLIMITED INDENTATION LENGTH, FIXED VOLUME 

Let us first consider the question: 

(iv) Given the wing, body length, and total volume 
of the combination, what is the area distribution of the 
body which yields a minimum wave drag if the apexes 
of the Mach forecone and aftercone enclosing the wing 
he within the body (see fig. 6) and the specified volume 
is large enough for the body to be real? 



multipole distributions (or the second body’s equivalent 
multipole distribution) in equation (6). Since only the 
sources interfere, the drag can be written in the form 


D=D, 


q C L ° / i 7 f'o , _ . . 

+n ■ (38 > 


where the interference term has been integrated by parts 
and — Lo', La, and — Z Q , la form bounds of the arbitrary and 
Sears-Haack source distributions, Oo, and a<) , respectively. 
As is well known 


and 


Dsa 8Vag 
q ~ vlo* 


(39) 




'(*)= 


8V sg 2x*-Z 0 2 

rio 4 


(40) 


Placing equation (40) in (38) and integrating, one finds 

Oo/Sl) 


D=D, 


8 g\W 2 


Vsu 

and since one can easily show 


o GoA* 


■ dx\-\-D\ 


v ^-mr~w Xl ^ Xl)dXl 


equation (37) follows immediately. 

'Returning now to problem (iv), we see that its solution 
follows from equation (37) and the solution is, in fact, simply 
a Sears-Haack body having, at the appropriate place relativo 
to the wing-body juncture, the additional area variation 
specified by equation (4) . This follows, since, if D\ represents 
the combined drag of the wing and indentation, then V it 
the combination volume of the wing and indentation, is zero. 
Hence, the minimum value of D, for a given volume, is 
obtained when D S h and D x are independently minimized. 
But Dss is already a minimum on a volume basis and D t is a 
minimum for a given wing. Notice the location of the 
wing along the body is immaterial, provided the required 
indentation can be accommodated by the fuselage. 


This problem can be solved in a simple manner by means 
of the following lemma discovered by R. T. Jones, see 
reference 1. 

Designate the closed body which has a minimum wave 
drag for a fixed volume and length as a Sears-Haack body. 
Then the total drag of a Sears-Haack body and any other 
wing or (slender) body entirely within the fore and after 
Mach cones with apexes at the tail and nose of the Sears- 
Haack body, respectively, is given by the equation 

D^D SB (l+^)+D 1 (37) 

where 

D slI wave drag of Sears Haaek body alone 

Di wave drag of second body alone 

Vsb volume of Sears-Haack body 
volume of second body 

A proof of this lemma can be obtained by placing the 
Sears-Haack source distribution and the wing equivalent 


LIMITED INDENTATION LENGTH ON 8EAB-HAACK BODY, FIXED VOLUME 

Consider, next, the more difficult problem 

(v) Given a wing and Sears-Haack body of length 
2la (long enough to contain the apexes of the fore and 
after Mach cones enclosing the wing), what modification 
of this fuselage within the length 1/ + (and within 
that length only, see fig. 7)- minimizes the total wave 
drag for a given total volume? 

In order to answer this question, it is necessary to consider 
separately two cases; namely, the one in which li>W and 
Zi>io (i.e., the portion of the body free for variation contains 
the apexes of the wing’s Mach cone envelope, as shown in 
fig. 7) and the other in which the preceding conditions are 
not satisfied. 

First consider the combination for which 1\>L </ and 
Li>Lq. The wave drag of such a combination can always 
be calculated using equation (37) wherein D su is the wave 
drag of tlie basic Sears-Haack body fixed by the stationary 
nose and tail portions, D t is the combined wave drag of the 
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Fioure 7 . — Extent of cancellation sources for wing and Sears- 
Haack body. 

Wing and the (os yet unknown) body indentation, and Vi is 
the net difference in volume between V SB , the volume of 
the Sears-Haack body, and the final volume of the complete 
configuration. Since the basic Sears-Haack body is fixed 
and the total volume is given, the entire term Dsb[1-\-(2ViI 
Fsw)] is fixed and the solution to the problem is obviously 
that for which the wing equivalent sources and the source 
simulating the body indentation combine to form a Sears- 
Haack distribution in the interval — li 
Using equations (2) and (3) to relate the wing and body 
source variations to their respective areas, we find the 
fuselage cross-sectional area can be written for — — W 

for —W<x<h 

■5^5, »+*)((.-*))>* (41b) 


and for l\<x<l Q 


S/x)= 


3 * 10 * 




(41c) 


The total wave drag of the wing and the fuselage, as given 
by equations (41) is then 


D 8V sb V, , , 12817 

2 ~ rio* V^VsJ^^+l/y 


De 


where De is defined by equations (5), (6), and (7). 

Since, as we have been assuming, PA is small, De is negli- 
gible, and a comparison between equations (37) and (42) 
shows that the drag of the combination formed by mounting 
two wing panels on a Sears-Haack body can be reduced 
without a change in the total volume and with a modifica- 
tion limited to the interval — Z/^arOi by the difference 
between the drag of the two panels flying alone and the 
drag of a Sears-Haack body having a length equal to Z/-f Z„ 
and a volume equal to the volume bf the two panels. So 
long as the points x= — Z/ and x=h do not lie off the basic 
body, and so long as the required indentation can be accom- 
modated, this result is independent of the wing’s fore-and-aft 
position. 

If the body modification is limited so that either Z/^Z,/ 
or k<U or both, the above results do not apply, since, 
in such cases, the second body — in the sense defined above — 
cannot be varied for x between — Z/ and — JV or D> and li or 
both, and its drag cannot, therefore, be reduced to that of 
an equivalent Sears-Haack body. The best modification in 
this case can be calculated from the results presented in the 
material immediately following. 

LIMITED INDENTATION ON ABpiTBABY BODY-FIXED VOLUME 

Consider the question 

(vi) Given a wing, a body length, and the area dis- 
tribution of the fore-and-aft portions of a body, what 
is the variation of area along the intermediate portion 
of the body which yields a minimum wave drag for a 
fixed total volume? 

Again, as in equation (17), let <r(z) represent the sum of the 



Figure 8. — Extent of modification for wing and basic body (not 

Sears-Haack). 


( 42 ) 
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sources representing the basic body and the wing equivalent 
source distribution, 


c{z)=s f (z)+± j* T s 0 (x,e)de 


(43) 


It is now convenient, however, to let <r(x) be a fixed function 
in the entire interval — U <x<U, see figure 8, and let the 
body modifications, which are to be optimized in the interval 
— Zj'CarCfi-be represented by AS/(x) which has the end con- 
ditions 

A^/(-Z 1 )=AS/(Z 1 )=0 

A J S r X-Z 1 )=ASXZ 1 )= 


0=0 -) 
)=0 J 


(44) 


The change in volume caused by the body modification, AV, 
is given by " 

AV=—J I xAS/(x)dx . ’ (45) 

The usual variational procedure leads directly to the 
integral equation 


fXo+XiZ,; — Zi<Xi<Zi (46) 


Pi AS/'ixojdz? __ Po <r''(xx)dx 

J-q x 2 — x* — J-4 x^—Xi 

where Xo and \ are fixed by the conditions given in equa- 
tions (44) and (45). Equation (46) is similar in form to 
equation (13) and its inversion can be obtained by use of 
methods similar to those for inverting the latter equation. 
Thus, the solution to equation (46) becomes for —li<x<k 

g{ J;-; 


S, 


/“ — [ 2/ s -Z 1 %-2(F+Ay)]+^'V 

*1 fo— x)-yJXi *— h* H h J 


where 


" J-V V^-Zi 2 Jh 


’ll Xi n a'(z 1 )dx l 


(47) 

(48) 


and V is the total volume of the wing and unmodified 
fuselage, that is 

y=-J J x</{x)dx (49) 

Equation (47) integrates to give 


ASXz)=-<KaO+^rs Hi x+l^JTT^- 


[2I t -l^I 0 -2(V+AV)W-^+E(x) 


(50) 


where 


H(x)= - f ~ h <f(x\) T tan" 1 ■ =M M = =+|~| dx 1 - 

- f* <r'(*i) ["tan -1 - -i 1 ... +^~| (&! (51) 

*■ J'. L V(ii*-®*)(»i*-ii*) 2 -J 

If D c is the drag of the original wing-body combination 
and Dag is the drag of a body of revolution having the 
same normal area distribution as the modification, then 


7? [W+*V)+h%-2I 2 ) 

2 2 2 Mi 

On the other hand 2? r# can be written 

y-T+3? {u M-w-w'+AW+i-} 


where if 


<?(*!,«»)= 


._L_ /iOZ 

*i— *a Vx^-Zi 2 


(52a) 


(52b) 


(63) 


D*/2 can be expressed as 

Y = 2~ [ J-;/ < ^ 1 J-V ^'(ZiV'fo^foiZa)— 

2 J" ( ‘ dx 1 J i ‘ v dx 2 <j'(x l )o"{xJ)G(x 1 ,x 2 )- 1 r 

J t ° cZXiJ^ (/x 2 cr / (Zi)ff''(Z 2 )G , (Zi ) Z 2 )J 


(64; 


LIMITED INDENTATION ON ARBITRARY BODY— FIXED DIAMETER 

As a final example in. this section, consider the questior 
(vii) Given a wing, a body length, and tho area dis 
tribution of the fore-and-aft portions of a body, wha 
is the intermediate variation of fuselage area that ha 
a given area at some intermediate station x—d, an( 
yields a minimum wave drag for the combination? 

Using the same definition for <r(x) as is given in equatioi 
| (43), and again designating the area modification as A S f (x) 
one can apply the same methods used to develop equation 
(21) and (46) and write the integral equation for AS f (x) z 
the form 


P> AS/(x^)dx 2 _ P» <A{x 2 )dx 2 
X 1 —X 2 J —h' Xi—Zi 


f Xoj — W<*\<d\ 


(6f 


where Xq and Xi are constants whose values depend upon til 
restraints. 

The solution to equation (56) can be written 

J i. x — 1 J l\\x— a 

( 6 ' 

and the three constants A, B, and are fixed by tho coi 
ditions: (1) continuous slope 


AS/(x)dx=0 


(2) the body area at x=h is unchanged 

AS/(x)dx=Q 


Jv 


and (3) the body area at d 2 is given 


Jv 


AS/(x)dx=ASf{di) 


(5£ 


(58 


(58 
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The final solution is 

A SJ.X) = - a (x)+H(x)-^- (x-drY cosh- ^^,-1- 

<rtf, )--&&)]} 

where 

_^l_ [a SMd+<tdd-H(dd-&^ (M+JoV)] 

(60) 

and /* and H(x) are defined in equations (48) and (51), 
respectively. 

The drag can be expressed either as 


(69) 


—=——^^+‘irOiASj{di) 
2 2 2 


(61a) 


where D„ is again the drag of the original unmodified com- 
bination and D^s is the drag of the modification alone, or as 

(Bch-C^JIT^d?) 3 (61b) 


where D* and C x are defined by equations (54) and (60), 
respectively, and B is given by 


B^j* On i&W-df-ZD) 


(62) 


Ames Aeronautical Laboratory 
National Advisory Committee for Aeronautics 
Moffett Field, Calif., Aug. 4, 1955 
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